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ABSTRACT 

Extensive optical and near-infrared (NIR) observations of the type lb supernova 
1999dn are presented, covering the first year after explosion. These new data turn this 
object, already considered a prototypical SNIb, into one of the best observed objects 
of its class. The light curve of SN 1999dn is mostly similar in shape to that of other 
SNelb but with a moderately faint peak (My = —17.2 mag). From the bolometric light 
curve and ejecta expansion velocities, we estimate that about 0.1 1M© of 56 Ni were 
produced during the explosion and that the total ejecta mass was 4 — 6M with a 
kinetic energy of at least 5 x 10 51 erg. The spectra of SN 1999dn at various epochs are 
similar to those of other Stripped Envelope (SE) SNe showing clear presence of H at 
early epochs. The high explosion energy and ejected mass, along with the small flux 
ratio [Call]/ [01] measured in the nebular spectrum, together with the lack of signa- 
tures of dust formation and the moderate metallicity environment is not inconsistent 
with a single massive progenitor (Mzams ^23 — 25 M©) for SN 1999dn. 

Key words: Supernovae and Supernova Remnants: general - Supernovae and Su- 
pernova Remnants: 1999dn 



(e.g. Hakobyan et al. 2008) and, in particular, by their asso- 
ciation to massive star formation regions (e.g. Van Dyk et 
al. 1999; Anderson & Jones 2008). An additional evidence 
of a link between the different subtypes of CCSNe came 
with the discovery of objects, called type lib, that meta- 
morphose from type II at early stages to type lb later on. 
The prototype of this subclass is SN 1993 J in M81, one of 
the best studied SNe at all wavelengths ever (e.g. Barbon et 
al. 1995; Richmond et al. 1996b). 

While early time spectra of CC-SNe can be very dif- 
ferent as a consequence of the different configurations of 
the progenitors at the moment of their explosions, late time 
spectra of all CC-SNe are consistently similar with strong 
emission lines of neutral and singly ionized O and Ca, in 
addition to H Balmer lines for SNII. The standard scenario 
is that SNIb descend from massive stars that have lost their 
H envelope through strong winds or mass transfer to a com- 
panion (Heger et al. 2003). If in addition to the H also the 
He envelope has been removed, then the SN will appear of 



1 INTRODUCTION 

In the Supernova (SN) taxonomy, type lb Supernovae 
(SNIb) are defined as the subclass of Core-Collapse explo- 
sions (CC-SNe) whose early-time spectra are characterized 
by strong Hel lines (e.g. Wheeler & Levreault 1985; Turatto, 
Benetti, & Pastorello 2007). CC-SNe are thought to descend 
from massive progenitors (M > 7 — 8 M©) and include also 
the most frequent, H-dominated type II SNe (SNII) as well 
as type Ic SNe (SNIc), which appear deprived of both H and 
He. The physical connection among these subtypes is pro- 
vided by their location almost exclusively in spiral galaxies 

* Based on observations collected at the European Organisa- 
tion for Astronomical Research in the Southern Hemisphere, 
Chile (ESO programmes 64.H-0604 and 65.H-0292), at the Ital- 
ian 3.58m Telescopio Nazionale Galileo and the William Herschel 
Telescope (La Palma, Spain), and at the Copernico telescope (Asi- 
ago, Italy). 

f E-mail: stefano.benetti@oapd.inaf.it 
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type Ic. For this reason SNIb, Ic and lib are often referred 
to as Stripped Envelope (SE) SNe. 

With the improved quality of data, the differences be- 
tween type lb and Ic have shown to be subtle, and classifi- 
cations often controversial. For instance, in SN 19941, early 
on considered as the prototypical SNIc, has been found pos- 
sible signature of He (Filippenko et al. 1995; Clocchiatti et 
al. 1996); the type lb SN 1999ex was characterized by weak 
optical Hel but strong Hel A10830, A20581 lines in the near- 
IR (Hamuy et al. 2002); the peculiar SN 2005bf (Folatelli et 
al. 2006) and SN 2008D (Mazzali et al. 2008; Modjaz et al. 
2009) underwent a metamorphosis from a broad line type 
Ic at early times to a typical type lb at later epochs. The 
sharp distinction between the two classes seems therefore to 
be replaced by a continuity of properties in He abundances 
and/or excitation. 

The study of stripped-envelope SNIb/c has received 
fresh impetus in the past decade because of the association 
of some of them, in particular the most energetic SNIc (hy- 
pernovae), with GRBs of long duration (Galama et al. 1998; 
Hjorth et al. 2003; Malesani et al. 2004). More recently a few 
SNIc have been associated with less energetic X-ray flashes 
(Fynbo et al. 2004; Modjaz et al. 2006; Pian et al. 2006). 
An X-ray flash was also detected in the type lb SN 2008D, 
which was attributed either to shock break-out at the star 
surface (Soderberg et al. 2008) or to the effect of mildly rel- 
ativistic jets due to the collapse of a 30 M© star to a black 
hole (Mazzali et al. 2008). 

Despite this renewed interest the objects with detailed 
observations remain few. In particular, it is not well ascer- 
tained if SNIb do exist as a distinct class or if there is a 
uniform distribution of objects from SNII to SNIc with de- 
creasing H (and He) content in the outer layers. In this con- 
text, SN 1999dn is interesting because it has been adopted 
several times in the past to describe the average properties 
of SNIb (e.g. Branch et al. 2002; Chornock et al. 2010). 

SN 1999dn was discovered by Qiu et al. (1999a) on two 
unfiltered CCD images obtained on Aug. 19.76 and 19.82 
UT, respectively, in the Wolf-Rayet galaxy NGC 7714. The 
SN position is R.A.(2000.0)=23^36 m 14 s .7; Dec(2000.0) = 
+02°09 / 08".8, 9.9"E, 9.4"S of the galaxy nucleus (Qiu et al. 
1999a), in a region of steep background variation (Fig. 1). 
The parent galaxy, NGC 7714, is classified as SBb peculiar, 
and identified by Weedman et al. (1981) as a prototypical 
starburst galaxy. A second SN (2007fo), has been recently 
discovered 2.5"W, 12.4"N of the galaxy nucleus (Khandrika 
& Li 2007) , which showed similarly prominent He I lines and 
was also classified as SNIb (Desroches et al. 2007). Another 
highly reddened SN candidate was detected 2 X/ W, 5 /; N of the 
galaxy nucleus on UKIRT archival K-band images taken on 
1998, Dec. 5 (but not in H band, i.e. (H-K)> 1.2; Mattila 
et al. 2002). 

A series of spectra of SN 1999dn has been taken soon 
after the discovery by different groups which classified the 
SN as type la (Qiu et al. 1999b) and as type Ic because of 
the lack or weakness of He lines (Ayani et al. 1999; Turatto 
et al. 1999). Few days later the Hel lines emerged and the 
SN was re-classified as a SNIb/c by Pastorello et al. (1999). 
A week later the Hel lines strengthened again making the 
spectra of this SN very similar to those of other type Ib's. 

Due to the early discovery, SN 1999dn was observed 
by several groups (Deng et al. 2000; Matheson et al. 2001). 
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Figure 1. SN 1999dn in NGC 7714 and the sequence of local 
reference stars (cfr. Tab. 1). The image is a V frame taken with 
TNG+Dolores on June 28, 2000. North is up, east to the left. 

It soon became one of the best-observed SNIb and a test 
case for extensive modeling (Deng et al. 2000; Branch et 
al. 2002; Ketchum, Baron, & Branch 2008; James & Baron 
2010). In this paper we present original data collected at La 
Silla, La Palma and Asiago, and analyze them together with 
published material. The joint set of observations gives good 
multicolor photometry and dense spectral sampling, starting 
one week before maximum up to over two months. The SN 
has been recovered at late time in imaging (Van Dyk et al. 
2003) and spectroscopy (Taubenberger et al. 2009). 



2 OBSERVATIONS 

UBVRIJHK' photometry of SN 1999dn was obtained at 
ESO-La Silla and TNG-La Palma. Ten photometric nights 
were used to calibrate a local sequence of stars through com- 
parison with photometric standard stars (Landolt 1992). 
In turn, the local sequence was used to calibrate observa- 
tions obtained during non-photometric nights. The magni- 
tudes of the local sequence, labeled in Fig. 1, are reported 
in Tab. 1. The estimated errors (mean standard deviation) 
are reported in parentheses. Due to the small field of view of 
Arnica only star 1 was always included in the NIR frames, 
which then is the reference to check the calibration of the 
JHK' photometry. The averages of the measurements of this 
stars in three nights are reported in the footnote of Tab. 
1 along with their standard deviations. The moderate dis- 
persions of the measurements and their consistency with 
the 2MASS Point Source Catalogue (difference < 0.2mag) 
(Skrutskie et al. 2006) add confidence to the photometric 
calibration. 

The new photometric measurements of the SN are listed 
in Tab. 2. Observations were obtained on 19 different epochs 
up to one year after explosion. Data reduction followed stan- 
dard procedures making use of a PSF fitting technique for 
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Table 1. Magnitudes of the stars of the local sequence identified 
in Fig. 1. The photometric errors are in parentheses (in units of 
1CT 2 mag). 



star U B V R I 



* J= 12.25 ± 0.05, H= 11.50 ± 0.10, K' = 11.53 ± 0.10 



Table 3. Spectroscopic observations of SN 1999dn 



Date 


phase* 


inst.** 


exp. 


range 


resol. 




(days) 




(min) 


(A) 


(A) 


25/08/99 


-2.3 


DF 


45 


3600-9000 


9 


28/08/99 


+0.8 


DF 


45 


3500-9000 


9 


03/09/99 


+6.8 


DF 


45 


3500-9000 


9 


09/09/99 


+12.5 


WHT 


45 


3200-7550 


3.5 


14/09/99 


+17.5 


EK 


30 


3470-7470 


18 


14/09/99 


+17.8 


EF2 


15 


3400-7450 


14 


3/11/99 


+67.6 


DF 


120 


3500-9000 


10 


4/11/99 


+68.6 


DF 


120 


3500-9000 


10 


31/08/00 


+370.7 


EF2 


60 


6000-10250 


12 



* relative to the estimated epoch of B maximum (JD=2451418). 
**See note to Table 2 for coding, plus WHT = WHT 4.2m tele- 
scope + ISIS and EK = Asiago 1.8m telescope + AFOSC. 



1* 


15.20(05) 


14 


2 


15.76(06) 


15 


3 


17.91(13) 


17 


4 


19.36(16) 


18 


5 


22.82: 


21 


6 


16.19(02) 


16 


7 


15.68(05) 


15 


8 


19.83(09) 


18 


9 


22.71(09) 


21 



.73(03) 13.84(01) 

.45(04) 14.65(02) 

.35(06) 16.38(02) 

.58(08) 17.62(03) 

.57(09) 19.95(07) 

.20(04) 15.44(05) 

.69(03) 15.03(02) 

.43(09) 17.02(04) 

.31(09) 19.72(04) 



13.43(03) 12.86(02) 

14.20(02) 13.80(03) 

15.83(01) 15.36(04) 

17.07(09) 16.56(06) 

18.73(03) 17.44(04) 

14.97(03) 14.55(06) 

14.64(02) 14.27(02) 

16.13(04) 15.47(11) 

18.70(05) 17.58(12) 



the SN measurements. The mean photometric errors, esti- 
mated with artificial stars experiments, are given in paren- 
theses. 

The journal of the spectroscopic observations (Tab. 3) 
gives for each spectrum the date of observation (col.l), the 
phase relative to the adopted maximum (cfr. Sect. 3; col.2), 
the instrument (col. 3), the exposure time (col. 4), the wave- 
length range (col. 5) and the resolution derived from the aver- 
age FWHM of the night-sky lines (col. 6). In order to improve 
the signal-to-noise ratio, the average of the two spectra taken 
on 3 and 4 September 1999 is shown in Fig. 6, after check- 
ing that there was negligible evolution. In some spectra the 
telluric absorptions have not been removed because of the 
lack of suitable standard stars with sufficient signal-to-noise 
ratio. In other cases, it has been impossible to remove the 
contamination by the underlying HII region whose lines can 
be either under- (e.g. day +370) or over-subtracted (day 
+0.8). The absolute flux calibration was verified against the 
B, V and R photometry. In case of disagreement, the spectra 
were corrected by applying a constant factor. In fact the slit 
was normally aligned along the parallactic angle, in order 
to minimize atmospheric differential light losses (Filippenko 
1982). 



3 PHOTOMETRY 

The UBVRIJHK' light curves of SN 1999dn are shown in 
Fig. 2. The figure includes also the R-band data from Fig. 27 
of Matheson et al. (2001) (shifted by 15.85 mag to match 
our observations) , the late time observations of Van Dyk et 
al. (2003) and the pre-discovery limit by Qiu et al. (1999a). 
Only in the R band the light curve has been well monitored 
before and around maximum. A low order polynomial fit of 
all points around the peak provides JD^ acc = 2451420.5 =b 
0.5 (30 Aug. 1999) at R max = 15.85 ± 0.05 mag, in good 
agreement with the estimate by Matheson et al. (2001) (31 
Aug. 1999). In the other bands the observations started a few 
days later and the uncertainties on the epochs of maxima are 
larger (cfr. Tab. 4). Following common convention we adopt 
the epoch of B maximum (JD 2451418.0) as reference which 
is 3.5 days before the reference epoch used by Matheson 
et al. (2001) and, later, by Branch et al. (2002), Ketchum, 
Baron, & Branch (2008) and James & Baron (2010). 

The peak of the light curve is asymmetric with a steep 



rise to maximum followed by a slow decline for about 10 
days and a faster decrease afterwards. The contrast between 
maximum and inflection point as well as the width of the 
maximum light seems to be color dependent, with shorter 
wavelengths having narrower peaks and larger magnitude 
differences. A color dependence is visible also in the decline 
rate during the late radioactive tail, longer wavelengths hav- 
ing steeper slopes (cfr. Tab. 4). Only three epochs of NIR 
photometry are available: two during the early post peak 
decline and one in the radioactive tail. Though scanty the 
NIR photometry has been very useful for the construction 
of the bolometric flux (Sect. 3. 2). 

A comparison of the colour evolutions of SN 1999dn 
with several other SNIb/c is illustrated in Fig. 3. All objects 
have been dereddened according to values reported in Tab 5. 
An overall similarity is found among all objects of this var- 
iegated class but in the pre-maximum phase. At maximum 
light the B-V colors are between +0.1 to +0.6 mag and be- 
come redder afterward, reaching a maximum value (0.9 to 
1.2 mag) at about 20d because of the cooling due to ex- 
pansion. Thereafter the B-V colors slowly turn bluer due 
to the progressive increase of the emission line strengths. A 
good coverage before maximum is available only for the he- 
lium rich SNe 1993J, 1999ex, 2008D and 2008ax, which show 
very different behaviors. Whereas the blue-band light curves 
of the type lib SN 1993 J showed an early peak attributed to 
the emergence of the shock breakout, SN 2008ax, the other 
SNIIb of the sample, showed evidence of its emergence only 
in the very early UV light curve (Roming et al. 2009), but 
not in the optical (Pastorello et al. 2008; Taubenberger et 
al. 2010) and thus the color curve had an opposite behav- 
ior, starting red and turning bluer. Already one week be- 
fore maximum the B-V curves of both objects reach similar 
values (B-V= +0.15 mag) and evolve similarly afterwards. 
The pre-maximum difference may be attributed to the fact 
that SN 2008ax had a more compact progenitor and a lower 
density wind than SN 1993 J (Chevalier & Soderberg 2010). 
An early peak in the blue light curves is also visible in SN 
2008D and possibly in SN 1999ex. Indeed both objects show 
early B-V colors comparable to SN 1993J. However, the blue 
color peak is very short in SN 1999ex and the evolution ex- 
tremely fast, with the B-V color reddening by over 1 mag 
in only 3 days and resembling the general trend described 
for SN 2008ax thereafter. The behavior of SN 1999ex leaves 
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Table 2. Photometric measurements for SN 1999dn. The photometric errors are in parentheses (in units of 10 2 mag). 



date 


JD* 


U 


B 


V 


R 


I 


J 


H 


K' 


instr. 


OK /no /on 

zo/uo/ yy 


K 1 A 1 C "71 
01410. / 1 


lo.yo(Ub ) 


lb.oU(Uo ) 


1 1 Q ( C\A \ 

lb. lo(U4 ) 


i c n<5 1 c\a \ 
10.yb(^U4 J 


lo.oz(Ub J 








JJr 


OQ /HQ /on 

zo/uo/yy 


ki /i i q ?n 
0141o. < y 


iD.yo^UD J 


lb.ol(Uo J 


lb. 1U(U4 ) 


10.oU(U4 J 


10. / z(Ub ) 








JJr 


no /no /on 

Uo/uy/yy 


014z4. / / 


1 / .zo(Ub ) 


1 ( .1)1(1)0 ) 


~\ R OA ( C\A \ 

lb.z4(U4 ) 


i c no /n/i "\ 
10.yz(^U4 J 


lo.bo(Ub J 








JJr 


n^j /no /no 

Uo/uy/yy 


51427.77 


17.bz(Ubj 


17.15(U5J 


i <^ on/ nc^ 
lb.o(J(U5) 


lo.yo(Ubj 


15.bb(Ubj 








JJr 


1 1 /no /on 

ii/uy/ yy 


ki /i qq /in 
014oo.4y 


lo.oy^ub j 


1 7 71 /flK^ 

1 / . / 1(1)0 ) 


ib.oy(^u4 j 


lb.ll(Uo J 


10.by(^U0 ) 










i /i /no /on 

i4/uy/yy 


514cJ5.7o 


lo.71(Ub) 


i 7 on /no^ 
17.oU(Ud J 


lb.b4(UcJ) 


lb.l5(UcJ J 


15.ol(U4) 








ATP 


15/09/99 


51436.50 


18.97(09) 


17.87(05) 


16.76(03) 


16.28(03) 


15.85(03) 








EF2 


18/09/99 


51439.77 


19.28(09) 


18.14(05) 


16.91(03) 


16.41(03) 


15.95(03) 








EF2 


20/09/99 


51441.58 












15.38(15) 


15.16(10) 


14.92(30) 


ARN 


03/10/99 


51455.48 












15.78(10) 


15.50(10) 


15.31(20) 


ARN 


07/10/99 


51459.47 


19.99(09) 


19.01(06) 


17.81(05) 


17.09(05) 


16.63(05) 








OIG 


02/11/99 


51485.35 


20.14(15) 


19.59(06) 


18.32(06) 


17.56(06) 


16.88(06) 








OIG 


03/11/99 


51485.58 




19.53(05) 


18.32(03) 


17.57(05) 










DF 


04/11/99 


51486.52 




19.47(06) 


18.30(03) 


17.58(04) 


16.95(05) 








DF 


10/12/99 


51523.33 


20.28(10) 


19.89(08) 


18.82(07) 


18.02(07) 


17.30(07) 








OIG 


27/12/99 


51540.32 




19.94(08) 


19.11(07) 


18.20(07) 


17.59(07) 








OIG 


27/12/99 


51540.37 












17.93(20) 


17.46(25) 


17.64(30) 


ARN 


28/06/00 


51723.69 






22.57(35) 


20.78(20) 










Dol 


31/08/00 


51788.60 






^ 23.15(70) 


21.93(50) 










EF2 



* 2400000+ 

DF = ESO/Danish 1.5m telescope + DFOSC 

OIG = Telescopio Nazionale Galileo + OIG CCD Camera 

EF2 = ESO 3.6m telescope + EFOSC2 

ARN = Telescopio Nazionale Galileo + ARNICA IR Camera 
Dol = Telescopio Nazionale Galileo + Dolores 

room for the presence of early, short shock breakout signa- 
tures also in other objects. 

The colors of SN 1999dn are available only from around 
maximum, but the spectral evolution in the week before (cfr. 
Sect. 4) indicates that the SED has become bluer and bluer 
from day -6 to maximin light, which suggests a color evolu- 
tion for SN 1999dn more similar to SNe 1999ex and 2008ax 
than to SN 1993J. After reaching a maximum value of B-V 
similar to all other SNIb/c, SN 1999dn remained unusually 
red up to over 100 days. The adoption of an higher extinc- 
tion correction might reduce the difference (cfr. Sect. 3.1). 

The R-I color evolutions of SE-SNe are smoother 
and remain confined between —0.3 mag at maximum and 
+0.7 mag later on. All objects show similar evolution and 
SN 1999dn is in the troop. The outlier here is SN 1993J, 
which remains bluer at all epochs. 

3.1 Reddening and absolute magnitudes 

The interstellar NalD lines arising both from our and the 
parent galaxy are well seen in the low-resolution spectra at 
about 5893 and 5945 A, respectively. Therefore the redden- 
ing of SN 1999dn cannot be neglected. Schlegel et al. (1998) 
give a galactic reddening Emw(B — V) = 0.052 mag in the 
direction to NGC 7714. We have measured the EWs of the 
NalD components in the spectra of SN 1999dn on the spec- 
trum of highest signal-to-noise and resolution (WHT on 9 
Sept., res. 3.5 A) and found an EW M w(NaID)~ 0.51 A and 
EWN77i4(NaID)~ 0.45 A. These values are in good agree- 
ment with the average of EW measurements performed on 
all other available spectra with no evidence of significant 
time evolution. Assuming that the gas/dust properties in 



Table 4. Main data of SN 1999dn 



position (2000. 0) a 
parent galaxy 
offset wrt nucleus 
adopted distance modulus 
SN heliocentric velocity 
adopted reddening 



23^36™ 14.81 s +02 o 09™08 s .4 
NGC 7714, SBb pec b with starburst c 
9.9"E 9.4"S 
fi = 32.95 ± 0.11 
2630 ± 150 kms -1 
Emw(B-V) = 0.052 d Etot(B-V) = 0.10 ± 0.05 



U 
B 
V 
R 
I 

uvoir 



peak time 
(JD-2451000) 
418 ± 2 

418 ± 2 

419 ± 1 
420.5 ± 0.5 e 

424 ± 1 
419.5 ± 0.5 



peak observed 
magnitude 
16.9 ± 0.1 
16.8 ± 0.1 
16.1 ± 0.1 
15.85 ± 0.05 
15.60 ± 0.05 



peak absolute 
magnitude 
-16.6 ± 0.5 
-16.6 ± 0.4 
-17.2 ± 0.4 
-17.35 ± 0.40 
-17.55 ± 0.35 
= 2.0 X 10 42 erg s~ 



rise to max 



- 12 days 



explosion day ~ 406 




16 Aug. 1999 






late time decline 


interval 


late time decline 


interval 




mag(lOOd) -1 


days 


mag(lOOd) -1 


days 


U 


0.37 


67 -105 






B 


0.81 


67-122 






V 


1.79 


67-306 


1.43 


67-122 


R 


1.44 


67-370 


1.17 


67-122 / 


I 


1.64 


67-515 


1.19 


67-122 


uvoir 


1.49 


67-122 






M( 56 Ni) 


0.11 


M 







M(ejecta) 
explosion energy 



4-6 M Q 



a- Qiu et al. (1999a), b- http://leda.univ-lyonl.fr, c- Weedman et al. 
(1981), d- Schlegel et al. (1998), e- 30 Aug. 1999, f- our data only 



NGC 7714 and the Galaxy are the same, we estimate a 
total redening E* t(B — V) = 0.098 mag. A similar value 
Et t(B — V) = 0.131 mag is obtained by using the average 
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Figure 2. UBVRIJHK' light curves of SN 1999dn. The left panel is a zoom around the epoch of maximum, while the whole evolution 
is shown in the right panel. The R band observations by Matheson et al. (2001) are reported as small filled circles, shifted to match our 
R magnitudes. The prediscovery unfiltered limit by Qiu et al. (1999a) is reported in the R band scale (open circle + arrow). Late time 
(> 500d) F550W and F814W observations by Van Dyk et al. (2003) are reported in the V and I scale, respectively, as large symbols. To 
guide the eye, the late time V and R observations have been fitted by straight lines, while a dashed line is drawn to show the slope of 
56 Co decay. 



relation between E(B-V) and EW(NalD) (Turatto, Benetti, 
& Cappellaro 2003). Finally, the match of the color curves 
of SN 1999dn to those other SE SNe (cfr. Sect. 3) suggests a 
consistent, though formally slightly higher reddening, E(B- 
V)~ 0.15. Throughout this paper we adopt for SN 1999dn a 
total reddening E tot (B - V) = 0.10 ± 0.05 mag (Ay = 0.31 
mag). 

Since there is no direct measurements of the distance 
to NGC 7714, we use the Hubble's law. From the wave- 
length of the interstellar NalD absorption features we derive 
a recession velocity of 2630 ± 150 kms~\ which is consis- 
tent with the heliocentric radial velocity reported by LEDA 1 
(2797 =b 16 kms -1 ). LEDA provides also a velocity corrected 
for the Local Group infall into the Virgo cluster vy%r — 2798 
kms -1 . Adopting Ho = 72 kms -1 Mpc -1 , we derive a dis- 
tance modulus fi = 32.95 mag which is used throughout this 
paper 2 . 



Adopting the above mentioned values for extinction and 
distance, we obtain the following absolute magnitudes at 
maximum M v = -16.6 ± 0.5, M B -16.6 ± 0.4, M v 

17.2 ± 0.4, M R = -17.35 ± 0.40 and Mi = -17.55 ± 0.35, 
where in the computation of the errors we adopted an un- 
certainty of ±250 kms -1 on the Hubble distance modulus 
due to possible peculiar motion. 

We note that, after rescaling to a common distance 
scale, our determination of My is 0.5 mag brighter than 
that derived by Richardson, Branch, & Baron (2006), but in 
good agreement with the mean value of Mr = —17.01 ±0.41 
given by Li et al. (2010) for a sample of six SNIb, SN 1999dn 
included. In comparison to the subsamples studied in Li et 
al., SN 1999dn remains about 0.80 mag brighter than their 
unweighted mean for the normal SE SNe (SNIc+Ib+IIb). 



1 http://leda.univ-lyonl.fr 

2 A similar value (/i = 32.84) is obtained adopting the distance 



relative to the Virgo cluster in the 220 model (Kraan-Korteweg 
1995) with D(Virgo)=15.3 Mpc (Freedman et al. 2001). 
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Table 5. SE SNe used for color and bolo metric luminosity comparison. 



SN Host Galaxy SN type Etot(B-V) (m-M) main references 



1999dn 


NGC 7714 


lb 


0.10 


32.95 


Sect. 3.1 of this paper 


19901 


NGC 4650A 


lb 


0.16 


32.90 


Elmhamdi et al. (2004) 


1993J 


NGC 3031 


lib 


0.30 


27.80 


Barbon et al. (1995); Richmond et al. (1996b) 


19941 


NGC 5194 


Ic 


0.45 


29.75 


Richmond et al. (1996a); Dessart et al. (2008) 


1999ex 


IC 5179 


lb 


0.28 


33.19 


Stritzinger et al. (2002) 


2004aw 


NGC 3997 


Ic 


0.37 


34.17 


Taubenberger et al. (2006) 


2005bf 


MCG +00-27-5 


Ib/c 


0.05 


34.46 


Folatelli et al. (2006) 


2007Y 


NGC 1187 


lb 


0.11 


31.13 


LEDA, Stritzinger et al. (2009) 


2007gr 


NGC 1058 


Ic 


0.09 


29.84 


Valenti et al. (2008a); Hunter et al. (2009) 


2008D 


NGC 2770 


lb 


0.66 


32.29 


Mazzali et al. (2008); Modjaz et al. (2009) 


2008ax 


NGC 4490 


lib 


0.40 


29.92 


Pastorello et al. (2008); Taubenberger et al. (2010) 
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Figure 4. UBVRIJHK' bolometric (filled circles) and UBVRI 
(open squares) light curves of SN 1999dn. The points on the ris- 
ing branch and those after 4 months past maximum are largely 
based on extrapolations (assuming constant colors) and should 
be regarded as uncertain. The slope of 56 Co to 56 Fe decay is also 
displayed for comparison. 



Figure 3. Top panel: Comparison of the B-V color curve of 
SN 1999dn with those of other SE-SNe. The curves have been 
dereddened with the values reported in Tab. 5 and a standard 
extinction law. Bottom panel: as above, but for the observed R-I 
color curve (r-i for SN 2007Y). 



3.2 Bolometric light curve 

With the available photometry it has been possible to build 
the UBVRIJHK' bolometric light curve up to lyr (Fig. 4). 
When magnitudes in a bandpass were not available in a 
given night, the values were linearly interpolated. The light 
curves with a short temporal extension, e.g. those in the 
NIR, were extrapolated assuming constant colors. For this 
reason the bolometric flux during the rising branch and 
those at the latest epochs are based mainly on R band ob- 
servations and should be regarded as uncertain. The pho- 
tometry was corrected for extinction using Ry = 3.1. At 
the effective wavelengths of each filter monochromatic fluxes 
were computed; these were then integrated from the U to 
K ; bands using the trapezoid approximation, and converted 
to luminosity. The peak of the bolometric light curve is 



reached between the V and the R maxima on JD^cc — 
2451419.5 ± 0.5 at a luminosity L bo i = 2.0 x 10 42 erg s" 1 . 

The contribution of the NIR bands is substantial at all 
the epochs in which NIR data are available {Ljhk /L UVO i r = 
0.5, 0.5 and 0.4 at day 20, 40 and 121, respectively). This is 
more than a factor two larger than the value derived for SNIb 
2007Y (Stritzinger et al. 2009) and for SN 2008D (Modjaz 
et al. 2009), but similar to values derived for SNIIb 2008ax 
(Taubenberger et al. 2010). In the photospheric phase, the 
optical+NIR SED deduced from SN 1999dn photometry is 
consistent with a blackbody energy distribution with tem- 
peratures as derived from optical spectra (cfr. Sect. 4). The 
decline rate between 67d and 122d is 1.49 mag(lOOd) -1 , 
close to the average value of type la SNe. This value is sig- 
nificantly larger than the decline rate of 0.98 mag (100d) _1 
predicted if all the energy from the decay of 56 Co into 56 Fe 
was fully thermalized. No significant slope variation is ob- 
served up to day 370 (7 uvo ir = 1-65 [67d-305d]). 

In Fig. 5 we compare the bolometric light curve of 
SN 1999dn with those of other SE SNe. Unfortunately not 
all have coverage from optical to NIR wavelengths. The 
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comparison has been done both for the extended uvoir 
(UBVRIJHK', top panel) and the optical-only (UBVRI, bot- 
tom panel) bolometric curves. 

The rise to maximum can be very different. The type lib 
SN 1993J shows an early bright spike about 20d before a 
broader secondary maximum. For the other SNIb caught 
early-on, SN 1999ex (bottom panel), the bolometric light 
curve shows just an hint of the shock breakout. 
For the peculiar SN 2005bf a slow rise to a first maximum is 
observed followed about 25d later by a second brighter max- 
imum. The type lb SN 2007Y and the type Ic SNe 19941, 
2004aw and 2007gr were not detected sufficiently early af- 
ter the explosion and do not show any feature in the rising 
branch. The first reliable photometry of SN 1999dn has been 
obtained 6 days before B maximum (8.5d before R maxi- 
mum), while the available pre-discovery limit does not put 
stringent constraints on the date of explosion. The direct 
comparison with other SNIb/c suggests a rise time (relative 
to B maximum) similar to SN 2007gr, i.e. 11.5 ± 2.0 days 
(Hunter et al. 2009). 

The contribution of the JHK' bands to the bolometric 
flux of SN 1999dn stands out: while in the UBVRI domain 
it has almost the same luminosity as the type Ic SN 2007gr, 
in the uvoir it outpasses it by about 0.12 dex. 
The asymmetric peak of SN 1999dn, with a relatively fast 
rise and slow decline, reminds of the behavior of SN 2007gr, 
but with even slower decline. After maximum the bolometric 
light curve remains relatively broad and closely resembles in 
shape that of SN 2004aw, but remains about 0.3 dex fainter 
at all epochs. This simple comparison indicates that adopt- 
ing the Arnett model (Arnett 1982) at the early times, when 
the diffusion approximation is valid, the amount of 56 Ni 
synthesized in the explosion of SN 1999dn is between 0.10 
and 0.15 M (Miv*(2007gr)=0.076 M , Hunter et al. (2009); 
Miv;(2004aw)=0.30 M , Taubenberger et al. (2006)). Broad 
light curves are indicative of large diffusion time and hence 
of either a large ejected mass or a small kinetic energy. The 
comparison of the expansion velocity of SN 1999dn (cfr. Sec. 
4.1) with other SE SNe suggests that the long diffusion time 
is due to a large ejecta mass of SN 1999dn, much larger 
than that of the narrow-peak SN 19941 (cfr. Fig. 5). A more 
detailed discussion of the explosion parameters is given in 
Sect. 5. 



4 SPECTROSCOPY 

The entire spectroscopic evolution of SN 1999dn from about 
1 week before maximum to over lyr is illustrated in Fig. 6, 
thanks to the spectra listed in Tab. 3 and those published 
by Deng et al. (2000) and Matheson et al. (2001). The SN 
evolution is therefore well sampled at all crucial phases. We 
stress again that the phases adopted here, are relative to the 
date of B maximum (JD 2451418.0) and differ from those 
used in previous papers (cfr. Sect. 3). 

The first spectrum (— 6d, Aug. 21) has been subject 
to deep scrutiny. Deng et al. (2000) have modeled it with 
SYNOW using a black-body temperature Tbb =7500 K and 
photospheric velocity v p h=16000 kins" 1 . Besides the broad 
absorptions at 4750 and 4950 A (restframe) due to Fell, and 
at about 8200 A due to Call IR, the spectrum is dominated 
by the Hel A5876 and A7065 which contribute to the strong 
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Figure 5. Comparison of the bolometric light curve of SN 1999dn 
with those of other type Ib/c SNe. The phase is from B maximum. 
The comparison based on the whole optical-NIR (uvoir) domain 
is in the top panel, that based on the UBVRI bands (uBgVri for 
SN 2005bf, UuBgVri for SN 2007Y) in the bottom panel. The 
SN distances and reddening are reported in Tab. 5. Note that 
with the adopted distance and reddening SN 19941 is among the 
brightest objects. Small misalignments in the epochs of maxima 
are due to different choices of the reference epochs. 



5620 A and 6820 A features. The broad absorption at about 
6200 A has attracted the attention of Deng et al. and other 
modelers and will be extensively discussed in Sect. 4.1. This 
spectrum has been modeled also by Branch et al. (2002) who 
used SYNOW with line optical depths varying as v" n and 
n — 8 instead of e~ v ^ Ve (v e = 1000 kins" 1 ) as in Deng et al. 
(2000). The spectrum was well reproduced with Tbb= 6500 
K and v p h= 14000km s _1 and the same contributing ions. 
Ketchum, Baron, & Branch (2008) modeled this and the 
other spectra around maximum with their non-LTE code 
PHOENIX trying both a standard solar and a three times 
higher metallicity. Although they adopted a smaller extinc- 
tion (E(B-V) =0.052 mag) and later explosion epoch (12 
days before Aug.31) with respect to the values that we now 
consider more realistic (cfr. Tab. 4), their fit is good also for 
the non-thermal Hel lines. The other features were identified 
with Call, Fell and OI. Similar results were obtained also by 
James & Baron (2010), who correct the spectrum only for 
galactic extinction (Baron, private communication). Their 
model parameters (T moc /=6000 K and vo = 11000km s _1 ) 
are slightly smaller than the values (Tbb= 6700 =b 200 K 
and v p h= 12800km s _1 ) we derive by fitting the de-reddened 
spectrum. Note that our Tbb is a color temperature, differ- 
ent from Tmod which is much closer to an effective temper- 
ature, while the afore mentioned SYNOW Tbb has little 
physical significance (Deng et al. 2000). 

Our new observations provide a good coverage of the 
epoch of maximum with spectra on days —2.3 and +0.8. 
At this phase the object has become bluer and hotter, reach- 
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Figure 6. The overall spectral evolution of SN 1999dn (including spectra by Deng et al. 2000 and Matheson et al. 2001 as thin lines). The 
wavelength is in the observer rest frame and the flux is not corrected for reddening. The ordinate refers to the top spectrum; other spectra 
are shifted downwards with respect to the previous by 1.2 x 10 -15 (second spectrum) and 0.6 x 10 -15 erg s _1 cm -2 A -1 (others). The 
spectra taken on Nov. 3 and 4 have been averaged, while for Sept. 14 only the EFOSC2 spectrum is displayed. Vertical lines correspond 
to the positions of Hel A4472, A4921, A5876, A6678, A7065 and A7281, with expansion velocities of 13000 kms" 1 at the earliest epoch 
and then slowing down with time. Each spectrum is labelled with the phase with respect to the B maximum (JD=2451418.0). For clarity 
the last three spectra have been multiplied by the factors reported in parentheses. Telluric features are marked with the Earth symbol. 



ing Tbb= 9100 and 7800 K, respectively. The line contrast is 
now higher with more pronounced absorption troughs and 
broad emissions. The Hel lines have become stronger so that 
the A6678 line is clearly identified in addition to the A5876 
and A7065 lines. In the blue strong H & K Call and Fell lines 
are visible which indicate photospheric expansion velocities 
Vph = 10400 and 10000 kms -1 , respectively. 

The next available spectrum on day +4.1, Aug. 31, 
though quite noisy, has been modeled using both SYNOW 
(T bb = 7000 K and v ph = 12000 kms" 1 , Deng et al. 2000; 
Vph = 10000 kms, Branch et al. 2002) and PHOENIX 
(Ketchum, Baron, & Branch 2008) giving an effective tem- 
perature of 6000-6500 K and v ph =10000 kms" 1 . The fea- 
tures are well fit by the same ions as on day -6 (Call, Fell, 
Hel and OI). The model with standard solar metallicity re- 
produces fairly well also the He lines both in absorption 
and in emission. Again the best fit temperature by James & 
Baron (2010) is lower (T mod =5250 K). Our black-body fit 
to the dereddened spectrum of day +4.1 provides Tbb°° 6400 
K which is lower than the values derived for the two brack- 
eting epochs (T bb = 7800 and 7000 K for day +0.8 and 6.8, 



respectively) possibly because of a poor instrument response 
calibration. The position of the Fe lines at about 5000 A cor- 
responds to v p h= 10000 kms" 1 , in good agreement with the 
velocities derived in the spectral modeling. 

Of much better quality (S/N> 100) is the spectrum ob- 
tained on day +6.8. The fit to the SED provides Tbb= 
7000K, and the expansion velocity from the Fell lines de- 
creases to v p h= 8600 kms" 1 . 

In the following week the evolution slows down. Two 
spectra are available on day +12.5 and +14, the former 
with an excellent signal-to-noise (S/N= 150) and resolu- 
tion (3.5 A), the latter with wider spectral range (Matheson 
et al. 2001). Our analysis provides a similar temperature, 
Tbb^ 5500 K in both spectra, but significantly different 
photospheric velocities (v p h= 6400 and 4900 kms" 1 , respec- 
tively). Both the Galactic NalD lines as well as those orig- 
inating in the host galaxy are well resolved, as mentioned 
in Sect. 3.1. Hel A5876 is very broad and possibly heavily 
contaminated by NalD. Also the lines at A6678 and A7065 
clearly stand out, whereas others (A4472 and A7281) are less 
pronounced (the former probably blended with Mgll). The 
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well developed Hel lines make the spectrum closely resemble 
that of prototypical Type lb supernovae. Ketchum, Baron, 
h Branch (2008) have studied the spectrum of day +14: 
again the solar metallicity model fits the observations bet- 
ter, though the Call and Hel absorptions are too strong. 
Till starts to contribute significantly to the Fell dominated 
region below 5000 A. 

Three spectra have been obtained on Sept. 14 (day 
+ 17.5 to +17.8) at different telescopes (cfr. Tab. 3 and 
Deng et al. 2000). The SEDs of our two spectra are simi- 
lar (T bb = 4800 K) but the S/N ratio of that obtained with 
EFOSC2 (plotted in Fig. 6) is definitely higher. The Deng 
et al. (2000) spectrum, significantly bluer (Tbb= 5700 K), 
has been used in the modeling both by Deng et al. (2000) 
and Ketchum, Baron, & Branch (2008). Again the standard 
composition model seems to be preferred though with too 
strong Hel and Call absorptions. The index n of the density 
profile p oc r -n was decreased from 13 to 10 to improve the 
fit. 

The subsequent spectra (day +21 and +42, i.e. 
Sept. 17 and Oct. 8) come from Matheson et al. (2001). 
At both epochs the main features and the SEDs have not 
changed significantly with respect to day +18. We measure 
Tbb~ 4900 K on the dereddened spectrum, in good agree- 
ment with the temperature adopted in the modeling (respec- 
tively 4800 and 4600 K, Branch et al. 2002; James & Baron 
2010). The measured expansion velocities are v p h= 5500 
and 4600 kms -1 , respectively. In the latter spectrum one 
can notice that the CaIR triplet starts to be resolved with 
an absorption at about 8360 A, corresponding to a velocity 
v(Call) ~ 7700 kms -1 . As the temperature of the SN has 
decreased to about 5500K, the region bluer than 5000 A is 
largely dominated by Till lines (Ketchum, Baron, & Branch 
2008). This explains why the high metallicity model fails to 
reproduce the observations causing too strong absorptions. 

The latest spectrum of the photospheric series (day 
+68) is the average of two spectra taken with the same 
instrument on two consecutive nights. The Hel lines are still 
visible as well as the underlying continuum (Tbb= 4400 K). 
Nevertheless, the progressive transition to the nebular phase 
can be recognized. In particular, the broad [OI] AA6300,6364 
and[CaII] AA7291-7323 ([Oil] AA7320,7330) features start to 
appear clearly. 

The characterizing features of the nebular spectra of 
SNIb/c are, indeed, [OI] and [Call] emission lines. These 
are well developed at the epoch of our last spectroscopic 
observation (day 371). Once deblended into its two compo- 
nents, the [OI] line has an overall gaussian profile centered 
slightly redward of the restframe position (6309A) with 
FWHM~ 4500 kms -1 , though the relatively poor signal-to- 
noise ratio still allows for the presence of additional struc- 
tures in the line profile (Taubenberger et al. 2009). In any 
case there is no evidence for profile distortion or a blueshift 
due to dust formation as detected in other objects (e.g. 
SN 19901, Elmhamdi et al. 2004) in the same [OI] line. 

4.1 The expansion velocities and the presence of 
Hydrogen 

Hydrogen lines at early epochs distinguish the spectra of 
SNIIb from other SE SNe. Actually, faint H lines have been 
revealed unambiguously in the type lb SN 2000H and, with 
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Figure 7. Zoom of the 4600 A (left) and 6200 A (right panel) spec- 
tral regions during the first weeks of evolution of SN 1999dn. The 
x axes are in expansion velocity coordinates with respect to the 
rest frame positions of H/3 and Ha, respectively. The phases rela- 
tive to B maximum light are indicated on the right axis. To guide 
the eye two dashed vertical lines are drawn in the spectra around 
maximum, corresponding to expansion velocities of about 16900 
and 12100 kms" 1 . After day +12.5 the broad 6200A feature is 
replaced by two narrower and weaker features which slightly slow 
down with time. Corresponding lines are drawn also in the left 
panel relative to H/3. The dotted line in the right panel is drawn 
at the position of Hel A6678. 



lower confidence, in other SNIb (Branch et al. 2002; James & 
Baron 2010). These authors concluded that H is present with 
low optical depths in SNIb in general and that it is located 
in a detached shell with velocities as high as 11000-13000 
kms -1 . In this context, some attention has been paid in the 
previous spectral analyses of SN 1999dn (Deng et al. 2000; 
Branch et al. 2002; Ketchum, Baron, & Branch 2008; James 
& Baron 2010) to a broad feature present in the spectra 
around maximum at about 6200 A. 

The feature seen in the first epoch (day -6) was at- 
tributed by Deng et al. (2000) to Ha or CII A6580, arising in 
a detached layer with velocity v=19000 or 20000 kms -1 , re- 
spectively. Also Branch et al. (2002) interpreted the feature 
as due to Ha from a detached layer of H at 18000 kms -1 . 
The possible alternative identification as Sill, proposed by 
Woosley & Eastman (1997) for SN 1984L, results in an ex- 
pansion velocity (7300 kms -1 ) significantly lower than that 
of any other ion and is, therefore, considered unlikely. Simi- 
lar identifications were proposed for the spectra of day +4.1 
(Ha, Branch et al. 2002, possibly blended with CII, Deng 
et al. 2000) and day +17 and +21 (CII, Deng et al. 2000, 
or Fell, Branch et al. 2002). The spectra were revisited by 
Ketchum, Baron, & Branch (2008). They found that at early 
times, when the 6200 A feature is strong and broad, an uni- 
form atmosphere of three times solar metallicity, devoid of 
any H, could provide a plausible explanation, making the 
feature a blend of Fell lines and Sill A6355. At later times 
the feature splits into multiple, distinct, weaker features, 
and solar metallicity fits better. However, they admit that 
higher metallicity in the outer envelope (as seen in the early 
time spectra) is difficult to explain, also in the light of 
our direct determination of a somewhat sub-solar 
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Figure 8. Zoom of the Fell (left) and Call IR (right panel) spec- 
tral regions during the first months of evolution of SN 1999dn. In 
the left panel, the approximate positions of Fell A4924, A5018 and 
A5169 are marked. The velocity reported on the abscissa refers to 
the A5169 line and ranges from -14000 to -5500 kms -1 from day 
-6 to day 68 past B maximum. In the right panel two dashed lines 
are drawn to guide the eye for the main component (right) and 
the possible high velocity feature (left) of Call. 

metallicity for the SN environment (c.f.r. Sect 5). 

The simpler interpretation of the 6200 A feature, therefore, 
seems to be the presence of H. The new, specific analysis 
by James & Baron (2010) confirms this identification and 
could suggest a H mass of Mh ^ 10 -3 M© in an outer shell 
of solar composition above the He core. 

With this work we add new, high signal-to-noise spec- 
tra around maximum, and readdress the issue of the possi- 
ble presence of H. In Fig. 7 we have zoomed into the Ha 
(right panel) and H/3 (left panel) regions of the best avail- 
able spectra. To guide the eye in the right panel we have 
drawn two vertical dashed lines corresponding to velocities 
of —16900 and —12100 kms -1 at the earliest epoch, assum- 
ing the Ha identification. The dotted line at lower velocity 
indicates the position of Hel A6678 (the same as in Fig. 6). 
The line is tilted to roughly match the velocity decrement. 
Contrary to Hel, the 6200 A feature appears at constant ve- 
locity (—16900 kms -1 ) up to maximum light, while a notch 
at about —12000 kms -1 might be present. Should these fea- 
tures be due to single lines and not to the conspiring effect 
of line blending (e.g. Fell and Sill), the constancy in expan- 
sion velocity is an indication of detached layers. Note that 
the blue wing of the strong Hel A5876 line extends outward 
to about —18000 kms -1 on day -6 and —15500 kms -1 on 
day +0.8, close to the estimated velocity of the fastest H 
layer. 

By day +6.8 the broad 6200 A feature has disappeared. 
Starting on day +12.5 down to day +21 two distinct, re- 
solved features spaced by about the same amount as before 
are clearly visible. They were noted by Ketchum, Baron, 
& Branch (2008) who suggested the same possible inter- 
pretation with Fell, Sill, CII and Till. The slower compo- 
nent, if identified with Ha, shows marginally higher velocity 
than Hel A6678 and slows down with time at the same rate 
(~ 1000 kms -1 in lOd). 

To check whether this pair of lines are due to H, we have 



drawn two lines at corresponding velocities on the left panel 
of Fig. 7 relative to H/3. Although the H/3 optical depth is 
expected to be significantly smaller than that of Ha, we note 
that possible signatures of H/3 are recognizable at the ex- 
pected positions in the spectra of days -2.3 and +12.5, the 
latter having the best signal-to-noise and spectral resolu- 
tion. Though not compelling, we consider this an additional 
evidence that H is present in the spectrum of SN 1999dn 
both in a detached layer above the photosphere and mixed 
with Hel. 

A 6250A feature was identified as Ha in the early 
spectra of the peculiar SN 2005bf at comparable velocity 
(rsj —15000 kms -1 ) on the basis of the detection of similar 
tiny features at the corresponding position of H/3 (Folatelli et 
al. 2006). Also Call and Fell lines had components with sim- 
ilar expansion velocity. For this reason we have investigated 
the possible presence of high velocity features for both these 
ions also in SN 1999dn. In Figure 8 the spectral regions of 
interest are plotted in analogy to those in Fig. 7. Broad and 
shallow Fell lines seem present early on also in SN 1999dn at 
about —14000 kms -1 , a velocity exceeding that measured 
for Hel (—12700 kms -1 ) and expected for the photosphere 
on day -6. The Fell lines are easily detected at smaller ve- 
locity in the following spectra. Strong, broad Call IR triplet 
absorption is detected at a velocity comparable to that of 
Hel. In the spectra of highest signal-to-noise ratio (day +6.8 
and +14) one may also see another component about 5000 
kms -1 faster but the first detection of high velocity features 
of Call starting at this epoch makes the identification un- 
likely. We conclude that, in addition to Ha, there is evidence 
of weak high-velocity Fell features before maximum. 

In Fig. 9 (top panel) we summarize the evolution of the 
line velocities as a function of time. The standard devia- 
tion of the velocities derived for Hel lines is a = 400 — 700 
kms -1 , possibly due to blending of various lines (e.g. NalD 
for A5876, contamination by telluric features for A7065, or 
by narrow Ha emission of the parent galaxy for A6678) . We 
estimate uncertainties of the same order for the possible H 
lines because of their weakness. The photospheric velocity 
of SN 1999dn is in agreement with those of other normal 
SNIb (Branch et al. 2002) and is well fitted by the power- 
law, v ph oc t -2/(n-1) with n = 3.6. As for other SNIb, Hel 
seems to be undetached at the first epoch but detached af- 
terward. The lowest velocity of the Hel layer is measured at 
about 6000 kms -1 . The faster H is detached at all epochs 
with velocity ranging between 17000 — 13000 kms -1 . The 
slower H component remains about 1500 kms -1 faster than 
Hel at all epochs. The two H components somehow bracket 
the H velocities of other SNIb reported in Fig. 23 of Branch 
et al. (2002). 

The bottom panel of Fig. 9 compares the Fell velocity 
of SN 1999dn with those of other SE SNe. SN 1999dn has 
an expansion velocity smaller than SN 1999ex, but higher 
than other objects, e.g. SN 2007Y. Only at early phases the 
velocity seems to deviate from the power-law fit by Branch 
et al. (2002). At these early phases the expansion veloc- 
ity is, in fact, more similar to that shown by the energetic 
SN 2008D. The overall normal velocity behavior indicates 
that the broad (slowly evolving) light curve of SN 1999dn 
after the peak is probably not due to a low expansion veloc- 
ity but to a larger ejected mass. 
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Figure 9. Top panel: velocity evolution of the Hel, Fell , Call IR 
and Ha lines of SN 1999dn in the photospheric phase. The crosses 
are the averages of the velocities derived from the minima of Fell 
lines, adopted as indicators of the photospheric velocity; the cir- 
cles the average velocities of Hel lines; the triangles (starred tri- 
angles) the velocity of the fast (slow) component of Ha, and the 
squares the velocity of the Call IR triplet. Bottom panel: com- 
parison of the velocity evolution of the Fell lines in several SE 
SNe. In both panels the curve is the power-law fit from Branch 
et al. (2002). 



4.2 SN 1999dn and other Stripped Envelope SNe 

We have compared the spectra of SN 1999dn with those of 
other SNe by means of GELATO, the spectra comparison 
tool developed by Harutyunyan et al. (2008) which compares 
input spectra with those present in our archive. Not unex- 
pectedly, the best match is always with those of other SE 
SNe. In Fig. 10 we show the comparison with a number of 
well studied objects during the pre-maximum and maximum 
phase. At these epochs the objects that match SN 1999dn 
best are the extensively studied SNIIb SN 2008ax (Pastorello 
et al. 2008; Chornock et al. 2010; Taubenberger et al. 2010) 
and the energetic SN 2008D (Mazzali et al. 2008; Modjaz et 
al. 2009). 

For the first spectrum of SN 1999dn (-6d) the closest 
match is with SN 2008ax on days -9 and -8 (Pastorello 
et al. 2008; Taubenberger et al. 2010). The strong Ha, 
which unequivocally marks the presence of H in SN 2008ax 
at this epoch, corresponds in position to the 6200 A fea- 
ture of SN 1999dn. Note, however, that the spectra of 
SN 2008ax that best match those of SN 1999dn are at an 
earlier phase, probably because of a faster expansion veloc- 
ity of SN 1999dn. At this epoch, the numerical match with 
SN 2008D is poorer because its expansion velocity at such 
epoch is even higher (Mazzali et al. 2008). There is also a 



general resemblance to SN 2007Y, but the expansion veloc- 
ity of SN 1999dn is significantly larger. 

Also at maximum SN 1999dn best matches SN 2008D 
and SN 2008ax before their B maxima. The SED of 
SN 1999dn, however, is slightly redder. The 6200 A feature 
corresponds to the fast blue wing of the Ha absorption of 
SN 2008ax. 

The comparison after maximum is shown in Fig. 11. On 
day +12.5 (top) the spectra of all objects are very similar. 
While Ha in SN 2008ax is still very strong, the 6200 features 
of both SN 1999dn and SN 2007Y have dimmed and two 
notches are left where the broad absorption was before (cfr. 
Fig. 7). Again, due to smaller expansion velocity (cfr. Fig. 9), 
the absorptions in SN 2007Y are redder than in SN 1999dn. 

A late-time (371 d) red-grism spectrum has been ob- 
tained for SN 1999dn (Fig. 11 bottom). In both SN 2008ax 
and SN 2007Y, a broad shoulder on the red edge of [OI] 
AA6300, 6364 is clearly visible, which might be interpreted 
as a sign of interaction of a fast (v~ 10000 kms -1 ) shell 
with circumstellar material. However, the late time interac- 
tion scenario is not fully consistent and other interpretations 
have been proposed, though not proved (cfr. Taubenberger 
et al. 2010). At the same position only unresolved Ha, [Nil] 
and [SII] lines from an underlying HII region are discernible 
in SN 1999dn, a sign that there is not (yet) CSM inter- 
action. [OI] AA6300, 6364 and [Call] AA7191, 7323 ([Oil] 
AA7300, 7330) are clearly visible, with a [CaII]/[OI] flux ra- 
tio of 0.55 ±0.10. 



5 DISCUSSION 

SN 1999dn was extensively studied in the last decade but 
our new data allow a better determination of the main pa- 
rameters (cfr. Tab. 4). The light curve results to be 0.5 mag 
brighter (My= —17.2) than previously estimated, but still 
fainter than the average of SNIb/c (M v = -17.70, H =72 
kms -1 Mpc -1 , Richardson, Branch, & Baron 2006). The 
multicolor observations (Sect. 3) have shown a relatively 
strong NIR flux. This large NIR flux is intrinsic, since the 
EW(NalD), the broad band colors and the SED all sug- 
gest SN 1999dn to be only mildly reddened (E(B-V)= 0.1, 
cfr. Sect. 3.1). The NIR accounts for about 50% of the 
flux during the advance photospheric phase. The bolometric 
light curve reaches the same luminosity as SNe 2007gr and 
1999ex, but is brigther than SN 2007Y, and has a width 
similar to that of SN 2004aw (cfr. Fig. 5). It is known that 
the peak luminosity gives hints on the ejected 56 Ni mass, 
while the width carries information on the total mass of 
the eject a and kinetic energy. To extract in a self consistent 
manner information on the ejected mass (M ej ), the nickel 
mass, (M( 56 Ni)), and the kinetic energy (Ex), we have em- 
ployed a toy model of the bolometric light curve (Valenti 
et al. 2008b). Our procedure is based on a two-component 
analytical model to account for the photospheric and neb- 
ular phases. During the photospheric phase (t< 30 d past 
explosion) a homologous expansion of the eject a, spherical 
symmetry, and a concentration of the radioactive 56 Ni ex- 
clusively in the core are assumed (Arnett 1982). At late 
times (t ^ 60 d past explosion), the model includes the 
energy contribution from the 56 Ni- 56 Co- 56 Fe decay, follow- 
ing the works of Sutherland & Wheeler (1984) and Cappel- 
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Figure 10. Comparison of the pre-maximum (-1 week, top) 
and maximum- light (bottom) spectra of SN 1999dn with those 
of other SE SNe. The spectra are plotted in the SN restframe 
and dereddened according to Tab. 5. The phase relative to the 
B maximum of each object is reported close to the SN identifier. 
The spectra come from Mazzali et al. (2008) and Valenti et al. 
(2011) for SN 2008D, from Pastorello et al. (2008) and Tauben- 
berger et al. (2010) for SN 2008ax, from Stritzinger et al. (2009) 
for SN 2007Y, and from Deng et al. (2000) and this paper for 
SN 1999dn. 

laro et al. (1997). The incomplete trapping of 7-rays and 
positrons has been accounted for using the Clocchiatti & 
Wheeler (1997) prescriptions for both photospheric and 
nebular phases. The bolometric light curve model suggests a 
56 Ni mass M( 56 Ni)= 0.11 M . A comparison with the values 
from Table 2 of Hunter et al. (2009), obtained with similar 
methods, shows indeed that M( 56 Ni) is marginally larger 
than in SN 2007gr (0.076 M ), and SN 2007Y (0.06 M , 
Stritzinger et al. 2009), but smaller than 2004aw (0.3 M ). 
Using an average opacity of 0.06 g/cm (as in Valenti et al. 
2008b) and a scale velocity v S c=v P h = 10100 kms -1 , we 
obtain an ejected mass M e j — 4 — 6 M and kinetic energy 
Er = 5.0 — 7.5 x 10 51 erg. Surprisingly, but not totally unex- 
pected because of the broadness of bolometric light curve of 
SN 1999dn, these values are comparable to those obtained 
for the energetic SN 2008D associated to an X-ray flash 
(Mazzali et al. 2008). In fact our toy model makes use of 
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Figure 11. As in Fig. 10 but about 2 weeks (top) and 1 
year (bottom) post-maximum. The spectra come from Mod- 
jaz et al. (2009) for SN 2008D, Taubenberger et al. (2010) for 
SN 2008ax, Stritzinger et al. (2009) for SN 2007Y and this paper 
for SN 1999dn. 



the scale velocity of the explosion (v sc =v P h, which are simi- 
lar for SNe 1999dn and 2008D) without taking into account 
possible differences in the density profiles. 

In the case of SN 1999dn this may cause an over- 
estimate of the physical parameters, the real values being at 
the lower boundary of those ranges. A more precise deter- 
mination should rely on more detailed codes. A comparison 
between the physical parameters deduced, for a sample of 
well studied SNIb, with our toy model with those derived 
by more sophisticated codes will be given in Valenti et al 
2010 (in preparation). Anyway, even considering the possi- 
bility that some physical values for SN 1999dn are slightly 
over-estimate, the supernova seems to belong to the group 
of relatively massive and energetic SNe lb. 

We noticed in Sect. 4.1 that the evolution of the photo- 
spheric velocity of SN 1999dn (as well as those of Hel and 
H) are in general very similar to those of other SNIb which 
implies, following Branch et al. (2002), comparable density 
profiles, masses and kinetic energies above the photospheres. 
As noticed by those authors, the strong similarity does not 
leave much room for asymmetries in SNIb, which is con- 
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firmed by the [01] line profile of SN 1999dn at late time 
(Taubenberger et al. 2009). Only the earliest (-6d) mea- 
surement of the Fell velocity is significantly larger than the 
power law fit by Branch et al.. 

The previously available spectra of SN 1999dn before 
and after maximum light have been extensively modeled by 
means of both the highly parametrized SYNOW (Deng et al. 
2000; Branch et al. 2002) and the more sophisticated non- 
LTE PHOENIX codes (Ketchum, Baron, & Branch 2008; 
James & Baron 2010). The latter successfully reproduces 
all main features of the spectra with a homogeneous stellar 
atmosphere of H, He, C, N, Ne, Na, Mg, Si, Ca and Fe 
with the main features identified as Hel, OI, Call and Fell 
at all epochs. Also the Hel absorptions were successfully 
reproduced including 7— ray deposited by the radioactive f3 
decay of 56 Co. 

Our new spectra confirm the presence of evolving fea- 
tures around 6200 A during the first weeks (Sect. 4). Up 
to maximum light the region is dominated by a broad and 
strong feature which, if attributed to H, implies expansion 
velocities of about 17000 krns" 1 (cfr. Sect. 4.1 and Fig. 7). 
At subsequent epochs (days +6.8 to +21) our new, supe- 
rior quality spectra confirm the presence of minor features 
(cfr. Ketchum, Baron, & Branch 2008) in the same region. 
Whether these components are due to H is not clear but the 
detection in the highest signal-to-noise spectra of weak ab- 
sorptions at the same velocities with respect to H/3 seems to 
support the existence of two H layers, one detached at veloc- 
ities v = 17000 to 13000 krns" 1 , comparable to that of H in 
other SNIb (Branch et al. 2002), the other only marginally 
faster than Hel (A(v) ~ 1500 krns -1 , cfr. Fig. 9) and possi- 
bly located in the outer He layer. The high Fell velocity at 
the earliest epoch (-6d) seems compatible with the existence 
of high-velocity H-rich layers. 

SN 2007Y and other SNIIb show evidence of interaction 
of the fast expanding H-rich layer with the CSM (Fig. 11, 
bottom) about one year past maximum light. This is not 
the case in SN 1999dn, indicating that no major mass loss 
episodes have occurred in the last decade before the explo- 
sion (assuming wind velocities ~ 2000 krns -1 , typical of 
WR progenitors). There is also no evidence of dust forma- 
tion at any time, neither from photometric (light and color 
curves) nor from spectroscopic observations (line positions 
and shapes). 

On the spectrum obtained at WHT on Sept. 9, 1999, 
we have measured the 03N2 index (Pettini & Pagel 2004) 
of the region adjacent to SN 1999dn along the slit of the 
spectrograph, corresponding to projected linear distance of 
~ 200 pc at the adopted distance of the SN. Relation (3) 
of Pettini & Pagel (2004) thus provides an oxygen abun- 
dance at the SN location 12 + log(0/H) = 8.39 + 0.05 + 0.25 
(where the first error is statistical and the second one is 
is the 95% spread of the calibrating relation), which is 
in excellent agreement with the estimate of Modjaz et al. 
(2010) (8.32) and slightly lower than the solar abundance 
12 + log(0/H) = 8.69 + 0.05 (Asplund et al. 2009) and close 
to the average metallicity derived for the sites of a sample 
of fifteen SNIb, which includes also SN 1999dn, studied in 
Modjaz et al. 

A search for the progenitor of SN 1999dn has been car- 
ried out on HST archive images (Van Dyk et al. 2003). The 
progenitor was not detected down to ^ — 7.3 mag and 



(U-V)o ^2.5 mag (these values do not change significantly 
with our assumptions on H and E(B-V)). Unfortunately 
this determination does not strongly constraints the nature 
of the progenitor either as a massive, single WR star (e.g. 
Mzams ^ 23 - 25 M at Z=0.02, Georgy et al. 2009), or as 
a star of lower initial mass in an interacting binary system 
(e.g. Mzams ^14-16 M at Z=0.02, Yoon, Woosley, & 
Langer 2010). However the lack of signatures of dust favors 
the single, massive star scenario, given the fact that while 
the radiation field of single WR stars is expected to prevent 
dust formation in their local environment, while binarity in 
WR stars seems to provide the necessary physical conditions 
for it (Crowther 2007, and references therein). 

The moderate metallicity environment (12 + 
log(0/H) ~ 8.4, slightly sub-solar) in which SN 1999dn 
exploded is not inconsistent with the single scenario. 
In fact, the binary channel for producing WR stars in 
WR galaxies (as NGC 7714) is important just at lower 
(12 + log(0/H) < 8.2) metallicity (e.g. Lopez-Sanchez & 
Esteban 2010), while the probability of forming single WR 
stars increases with the metallicity because it is easier to 
reach the WR phase due to the metallicity-dependence 
of the stellar wind (e.g. Lopez-Sanchez & Esteban 2010, 
and references therein). The lower mass limit for having a 
WR progenitor decreases from ^ 40 — 50 M at Z=0.004 
to - 23 - 25 M at Z=0.02 (Georgy et al. 2009). As a 
consequence the progenitor of SN 1999dn could be a single 
WR star having a main sequence mass ^ 23 — 25 M . 

Also the relatively small flux ratio [CaII]/[OI]= 0.55 ± 
0.10, known for being constant with time (cfr. the spectra 
collection by Taubenberger et al. 2009) is consistent with a 
single massive star scenario. In fact, Fransson & Chevalier 
(1987, 1989) have shown that this ratio is a diagnostic of the 
core mass of the progenitor, with higher ratios indicative of 
smaller cores. The ratios measured for SN 1993J, SN 2007Y 
and SN 2008ax are 0.6, 1.0 and 0.9, respectively (Stritzinger 
et al. 2009; Taubenberger et al. 2010). Since the core mass 
is strongly dependent on the progenitor ZAMS mass, thus 
we have an indication that the progenitor of SN 1999dn is 
more massive than the above mentioned SNe. 



6 CONCLUSIONS 

We have presented detailed photometric observations and 
new spectra of SN 1999dn from before maximum to the 
nebular phase. These new data turn this object, already con- 
sidered a prototypical SNIb, into one of the best observed 
objects of this class. 

SN 1999dn was a moderately faint SNIb (My = -17.2 
mag) which produced 0.11 M of 56 Ni. With a toy model 
we have estimated an ejected mass of 4-6 M with Er = 
5.0 — 7.5 x 10 51 erg. Due to the rough approximation of 
the model, these values may be slightly over-estimated. Our 
analysis on SN 1999dn confirms that, contrary to early be- 
lief, a prototypical SN lb may produce several foe of kinetic 
energy and eject several solar masses. 

Overall the main parameters of the explosion are com- 
parable to those of the type Ic SN 2004aw and the massive 
type lb SN 2008D, much higher than those of the low-energy 
and low-ejected-mass SN 2007Y. Higher explosion energy 
and ejected mass, along with the small flux ratio [CaII]/[OI], 
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the lack of signatures of dust formation and the relatively 
high-metallicity environment point toward a single massive 
progenitor (Mzams ^ 23 — 25 M ). On the other hand, none 
of these evidences completely rule out the scenario of a less 
massive star in a binary system. 

The spectra of SN 1999dn at various epochs are sim- 
ilar to those of other SE SNe that show clear presence of 
H at early (type lib SNe 1993J and 2008ax, type lb SNe 
2000H, 2007Y and 1999ex) or late (SNe 1993J, 2008ax, 
2007Y) epochs. Such similarities, coupled to the fact that ac- 
curate spectral modeling (e.g. Ketchum, Baron, Sz Branch 
2008) did not find other satisfactory explanations for the 
puzzling 6200A feature, lead us to support its identification 
with detached Ha. We conclude, therefore, that it is likely 
that residual H can be recognized in the spectra of most 
SNIb if observed sufficiently early on. 



ACKNOWLEDGMENTS 

The SN 1999dn spectra published by Deng et al. (2000) and 
Matheson et al. (2001) have been retrived from the Suspect 
database ( 

http : / /bruf ord . nhn . ou . edu/~ suspect/ index 1 . html 

) . We acknowledge the usage of the HyperLeda database ( 
http : //leda . univ-lyonl . f r 

) 

We thank the referee, Chornock, R., for his helpful com- 
ments. SB, MT, EC and FB are partially supported by the 
PRIN-INAF 2009 with the project "Supernovae Variety and 
Nucleosynthesis Yields" . 



REFERENCES 

Anderson, J.P., James, PA. 2008, MNRAS, 390, 1527 
Arnett W.D., 1982, ApJ, 253, 785 

Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, 

ARA&A, 47, 481 
Ayani, K., Furusho, R., Kawakita, H., Fujii, M., Yamaoka, 

H., 1999, IAUC 7244, 1 
Barbon R., Benetti S., Cappellaro E., Patat F., Turatto 

M., Iijima T., 1995, A&AS, 110, 513 
Branch D., et al., 2002, ApJ, 566, 1005 
Cappellaro E., Mazzali P. A., Benetti S., Danziger I. J., 

Turatto M., della Valle M., Patat F., 1997, A&A, 328, 

203 

Chevalier, R. A., & Soderberg, A. M. 2010, ApJL, 711, L40 

Chornock R., et al., 2010, arXiv: 1001. 2775 

Clocchiatti A., Wheeler J. C, Brotherton M. S., Cochran 

A. L., Wills D., Barker E. S., Turatto M., 1996, ApJ, 462, 

462 

Clocchiatti, A., Wheeler, J.C. 1997, ApJ, 491, 375 

Crowther P. A., 2007, ARA&A, 45, 177 

Deng, J.S., Qiu, Y.L., Hu, J.Y., Hatano, K., Branch, D., 

2000, ApJ 540, 452 
Desroches, L.B., Wang, X., Ganeshalingam, M., Filip- 

penko, A., 2007, CBET 1001 
Dessart L., et al., 2008, ApJ, 675, 644 



Elmhamdi A., Danziger I. J., Cappellaro E., Della Valle 
M., Gouiffes C, Phillips M. M., Turatto M., 2004, A&A, 
426, 963 

Filippenko A. V., 1982, PASP, 94, 71 

Filippenko A. V., et al., 1995, ApJ, 450, Lll 

Fransson C. & Chevalier R., 1989, ApJ, 343, 323 

Fransson C. & Chevalier R., 1987, ApJ, 322, 15 

Freedman W. L., et al., 2001, ApJ, 553, 47 

Folatelli C, et al., 2006, ApJ, 641, 1039 

Fynbo J.P.U., et al., 2004, ApJ, 609 962 

Galama T. J., et al., 1998, Nature, 395, 670 

Georgy C, Meynet G., Walder R., Folini D., Maeder A., 

2009, A&A, 502, 611 
Hakobyan A. A., Petrosian A. R., McLean B., Kunth D., 

Allen R. J., Turatto M., Barbon R., 2008, A&A, 488, 523 
Hamuy M., et al., 2002, AJ, 124, 417 
Harutyunyan A. H., et al., 2008, A&A, 488, 383 
Heger A., Fryer C. L., Woosley S. E., Langer N., Hartmann 

D. H., 2003, ApJ, 591, 288 
Hjorth J., et al., 2003, Nature, 423, 847 
Hunter, D. et al, 2009, A&A, 508, 371 
James, S. & Baron, E., 2010, ApJ, 718, 957 
Kandrika, H., Li, W., 2007, CBET 997 
Ketchum W., Baron E., Branch D., 2008, ApJ, 674, 371 
Kraan-Korteweg R. C, 1995, yCat, 7098, 
Landolt, A.U., 1992, AJ 104, 340 

Li, W. et al., 2010, MNRAS submitted, arXiv:1006.4612 
Lopez-Sanchez A. R., Esteban C, 2010, arXiv: 1004.0051 
Malesani D., et al., 2004, ApJ, 609, L5 
Matheson, T., Filippenko, A.V., Li, W., Leonard, D.C., 

2001, ApJ 121, 1648 
Mattila S., Meikle P., Walton N., Greimel R., Ryder S., 

Alard C, Lancon A., 2002, IAUC, 7865, 2 
Mazzali P. A., et al., 2008, Sci, 321, 1185 
Modjaz M., et al., 2006, ApJ, 645, L2 
Modjaz M., et al., 2009, ApJ, 702, 226 
Modjaz M., et al., 2010, ApJL, submitted, 

arXiv:1007.0661vl 
Pastorello, A., Turatto, M., Rizzi, L., Cappellaro, E., 

Benetti, S., Patat, F., 1999, IAUC 7245 
Pastorello A., et al., 2008, MNRAS, 389, 955 
Pettini M., Pagel B. E. J., 2004, MNRAS, 348, L59 
Pian E., et al., 2006, Nature, 442, 1011 
Qiu, Y.L., Qiao, Q.Y., Hu, J.Y., 1999, IAUC 7241 
Qiu, Y.L., Huang, 1., Yao, B., Li, H., 1999, IAUC 7244 
Roming, P. W. A., et al., 2009, ApJL, 704, L118 
Richardson D., Branch D., Baron E., 2006, AJ, 131, 2233 
Richmond M. W., et al., 1996, AJ, 111, 327 
Richmond M. W., Treffers R. R., Filippenko A. V., Paik 

Y, 1996, AJ, 112, 732 
Schlegel, D.J., Finkbeiner, D.P., Davis, M., 1998, ApJ 500, 

525 

Soderberg A. M., et al., 2008, Nature, 453, 469 
Stritzinger M., et al., 2002, AJ, 124, 2100 
Stritzinger M., et al., 2009, ApJ, 696, 713 
Skrutskie M. F., et al., 2006, AJ, 131, 1163 
Sutherland P. G., Wheeler J. C, 1984, ApJ, 280, 282 
Taubenberger, S., et al., 2006, MNRAS, 371, 1459 
Taubenberger, S., et al., 2009, MNRAS, 397, 677 
Taubenberger, S., et al., 2010, MNRAS, submitted 
Turatto, M., Rizzi, L., Salvo, M., Cappellaro, E., Benetti, 
S., Patat, F., 1999, IAUC 7244 



Supernova 1999 dn 

Turatto M., Benetti S., Cappellaro E., 2003, fthp.conf, 200 
Turatto M., Benetti S., Pastorello A., 2007, AIPC, 937, 187 
Valenti S., et al., 2008, ApJ, 673, L155 
Valenti S., et al., 2008, MNRAS, 383, 1485 
Valenti S., et al., in preparation 

Van Dyk S. D., Peng C. Y., Barth A. J., Filippenko A. V., 

1999, AJ, 118, 2331 
Van Dyk S. D., Li, W., Filippenko A. V., 2003, PASP, 115, 

1 

Weedman, D.W., Feldman, F.R., Balzano, V.A., et al., 

1981, ApJ 248, 105 
Wheeler J. C, Levreault R., 1985, ApJ, 294, L17 
Woosley S. E., Eastman R. G., 1997, thsu.conf, 821 
Yoon S.-C, Woosley S. E., Langer N., 2010, 

arXiv: 1004.0843 



